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Further  work  on  the  radiative  behavior  of  rocket  exhaust  plumes  at  high 
altitude  is  reported.  Comparisons  of  calculated  flow  fields  with  wind- 
tunnel  test  data  indicate  that  agreement  is  generally  satisfactory  if  allow¬ 
ances  are  made  for  absolute  calibrations  of  the  data.  Calculations  of  a  plume 
emission  according  to  some  proposed  mechanisms  have  been  performed.  It  appears 
that  production  of  OH  by  0  atom  attack  on  secondary  H  atoms  in  hydrocarbon 
fragments  can  result  in  substantial  chemiluminescent  emission  from  HjO  formed 
subsequent  reaction  of  the  OH,  Some  semi-quantitative  analogies  between 
rocket  exhaust  flows  and  material  releases  in  the  upper  atmosphere  are  dis¬ 
cussed,  and  derivations  directed  toward  comparing  collision  energy  distribu¬ 
tions  at  very  high  altitudes  are  presented. 
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ABSTRACT 

Further  work  on  the  radiative  behavior  of  rocket 
exhaust  plumes  at  high  altitude  is  reported.  Comparisons 
of  calculated  flow  fields  with  wind-tunnel  test  data 
indicate  that  agreement  is  generally  satisfactory  if  allow¬ 
ances  are  made  for  absolute  calibrations  of  the  data.  Cal¬ 
culations  of  plume  emission  according  to  some  proposed 
mechanisms  have  been  performed.  It  appears  that  production 
of  OH  by  0  atom  attack  on  secondary  H  atoms  in  hydrocarbon 
fragments  can  result  in  substantial  chemiluminescent  emission 
from  H2O  formed  by  subsequent  reaction  of  the  OH.  Some  semi- 
quantitative  analogies  between  rocket  exhaust  flows  and 
material  releases  in  the  upper  atmosphere  are  discussed,  and 
derivations  directed  toward  comparing  collision  energy  distri¬ 
butions  at  very  high  altitudes  are  presented. 
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This  research  was  supported  by  the  Defense  Advanced 
Research  Projects  Agency  of  the  Department  of  Defense  and 
was  monitored  by  Air  Force  Cambridge  Research  Laboratories 
under  Contract  No.  F19628 -72-C-0006 .  This  Second  Semi¬ 
annual  Technical  Report  coders  the  period  from  January  1, 
1972,  to  June  30,  1972. 

The  author  wishes  to  acknowledge  the  contributions 
of  Drs.  E.  R.  Fisher,  J.A.L.  Thomson,  and  I.  Oppenheim 
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FIGURE  CAPTIONS 


Vibrational  temperature  of  as  a  function  of  altitude 
in  the  ionosphere. 

Cumulative  radiant  intensity  of  the  "turbulent"  Apollo 
S-Ij  plume  at  120  km  in  the  H20  or  bands  upstream 
of  a  given  axial  distance  for  two  assumed  values  of  the 
N2  vibrational  temperature  according  to  the  "two-step" 

V-V,  V-T  mechanism. 

Cumulative  production  of  vibrationally  excited  N2  in  the 
"turbulent"  Apollo  S-II  plume  at  120  km  upstream  of  a 
given  axial  distance. 

Cumulative  radiant  intensity  of  the  laminar  Atlas  sus- 

tainer  plume  at  120  km  in  the  H20  v1  or  bands  upstream 

of  a  given  axial  distance  for  the  OH  chemiluminescence 

mechanism.  For  curve  I,  OH  is  formed  only  by  0  +  H2;  for 

curve  II,  the  reaction  0  +  C  H0  , 0  also  contributes.  Tur- 

n  zn+z 

bine  exhaust  composition  by  moles  is  26%  H2 ,  11%  CO,  4% 

C02,  1.3%  H2CO,  2.6  %  CH4,  2.6%  C^Hg,  19.5%  CnH2n+2. 

Cumulative  radiant  intensity  of  the  laminar  Atlas  sustainer 

plume  at  120  km  in  the  H20  or  v-j  bands  upstream  of  a 

given  axial  distance  for  the  OH  chemiluminescence  mechanism. 

Conditions  as  in  Figure  4  except  that  CH20  mole  fraction  is 

0.65%,  C  H0  . _  mole  fraction  is  20.2%. 
n  zn+z 

General  structure  of  rocket  exhaust  plumes  (a)  and  material 
releases  (b)  in  the  upper  atmosphere  at  low  ambient  pres¬ 
sures.  Continuum  tydrodynamics  assumed. 


7.  Coordinate  system  for  describing  rocket  exhaust  plume 
flow  and  interactions  with  the  free  stream. 

8.  Coordinate  system  for  describing  spherical  releases  and 
their  interactions  with  ambient  gas. 

9.  Surfaces  of  constant  first  collision  energy  between  at¬ 
mospheric  molecules  and  released  gases  when  the  initial 
velocity  of  the  released  material  is  large. 
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1 . 0  Introduction 

In  the  first  semi-annual  technical  report  on  this  con¬ 
tract  we  discussed  a  number  of  modelling  studies  of  high- 
altitude  (100  -  300  km)  exhaust  plume  emission,  predomi¬ 
nantly  in  the  near  and  middle  infrared,  and  concentrating 
on  processes  taking  place  in  the  mixing  region  between  plume 
gases  and  the  atmosphere.  It  is  generally  believed  that 
this  part  of  the  plume  is  responsible  for  the  emission  from 
large  vehicle  exhaust  plumes  measured  by  various  means  in 
the  near  infrared  region.  This  belief  has  been  confirmed 

for  at  least  one  vehicle  by  the  recent  ARPA-sponsored  AFCRL 

2 

me.  moments  of  the  Apollo  14  and  15  S-II  stages.  The  out¬ 
come  of  the  studies  described  in  our  earlier  report  was  that 
there  appeared  to  be  no  mechanism  or  set  of  mechanisms  yet 
investigated  which,  when  combined  with  rather  detailed  flow 
field  calculations  and  the  available  information  on  rates  or 
cross  sections  of  the  pertinent  reactions  and  excitation  pro¬ 
cesses,  produces  a  calculated  radiant  intensity  as  large  as 
that  observed  for  the  Atlas  and  Apollo  vehicles. 

During  the  period  covered  by  the  present  report,  we  have 
continued  these  modelling  studies.  The  activities  have  included 
a  fairly  detailed  examination  of  the  results  generated  in  the 
first  reporting  period,  a  comparison  of  the  fluid-dynamical 
part  of  the  calculation  with  experimental  data?  and  a  re-ex- 
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amination  of  some  of  the  chemistry  and  excitation  mechanisms 
included  in  the  models.  The  results  of  these  investigations 
can  be  summed  up  as  follows: 

1.  The  results  of  the  earlier  modelling  calcula¬ 
tions  appear  reasonable,  for  the  conditions 
examined. 

2.  The  fli id-dynamical  calculations  agree  reason¬ 
ably  well  with  experimental  measurements. 

3.  The  emission  produced  by  chemi-excitation  in  the 
exhaust  plumes  of  hydrocarbon  vehicles  is  sub¬ 
stantially  increased  when  0  atom  reactions  with 
unburned  hydrocarbons  are  considered,  and  this 
mechanism  could  be  responsible  for  a  substantial 
portion  of  the  radiant  intensity  of  the  Atlas 
plume- . 

A  number  of  additional  flow-field  calculations  have  been  per¬ 
formed;  as  yet,  most  of  these  have  not  been  investigated  for 
radiation  produced. 

In  addition  to  the  above  work  on  the  modelling  studies, 
work  ha.'  been  initiated  on  other  tasks  in  the  contract.  Re¬ 
visions  to  the  Exhaust  Plume  Phenomenology  Handbook?  prepared 
under  a  previous  contract,  have  begun.  We  anticipate  that  a 
revised  version  of  this  Handbook,  including  available  results 
from  the  ARPA  Plume  Physics  Program,  will  be  issued  later  this 
year.  Studies  of  simulation  of  exhaust  plumes  have  also  been 
initiated.  These  studies  are  predominantly  aimed  at  defining 
the  degree  of  simulation  of  exhaust  plume  processes  by  atmo¬ 
spheric  material  releases. 

These  studies  are  discussed  in  greater  detail  in  the  fol¬ 
lowing  sections. 


2 


«.ri*:*w.'Arv 


’W-JTK^T-'  T,p' 


2.0  Modelling  Studies 


2 . 1  Comparisons  with  experiment: 

As  a  means  of  validating  a  portion  of  the  modelling 

activity  conducted  under  the  first  half  of  the  contract,  we 

undertook  a  comparison  between  a  flow  field  experiment  con- 

3 

ducted  at  Arnold  Engineering  Development  Center  and  a  calcu¬ 
lation  conducted  here.  This  calculation  will  be  the  subject 

5 

of  a  separate  report,  so  only  the  high  points  will  be  inclu¬ 
ded  here. 

The  experiment  was  conducted  in  Tunnel  M  of  the  Von  Kar- 
man  Gas  Dynamics  Facility  and  involved  a  simulated  plume 
(represented  by  a  jet  of  heated  helium)  mixing  into  a  simu¬ 
lated  free  stream  (represented  by  a  surrounding  flow  of  N2) . 

The  forward  portion  of  the  plume  is  expected  to  be  in  the 

merged  layer  flow  regime.  Measured  data  include  dynamic 
2 

pressure  (pu  )  and  species  number  densities  (n^)  of  N2  and  He 
in  the  mixing  zone.  There  remains  some  question  as  to  the 
absolute  values  of  these  quantities  due  to  the  numerous  cor¬ 
rections  necessary  to  reduce  the  data.  In  addition,  rotational 
(Tr)  and  vibrational  (Ty)  temperatures  of  N2  were  measured, 
but  the  results  for  Ty  are  characteristic  of  reservoir  condi¬ 
tions  rather  than  those  in  the  test  section,  and  the  values 
of  vai*y  strongly  with  upper  3tate  J  level  and  permit  no 
easy  interpretation. 
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The  calculations  were  performed  with  the  MULTITUBE  and 
BOW  codes  for  the  specified  experimental  conditions.  The 
effect  of  the  forebody  housing  the  helium  nozzle  is  included, 
and  laminar  mixing  of  He  and  N2  is  assumed.  Comparisons  of 
the  calculation  with  the  experiment  indicates  that: 

2 

1.  Good  agreement  with  measured  values  of  pu  is  ob¬ 
tained  except  within  the  shocks.  (The  codes  assume 
the  shocks  to  be  discontinuities.) 

2.  Good  agreement  with  the  shape  of  the  measured  dis¬ 
tributions  of  nHe  is  obtained  (again  excepting 
values  at  the  shocks) ,  but  the  magnitudes  of  the 
measured  values  tend  to  be  40-60%  higher  than  the 
calculations.  This  discrepancy  is  thought  to  have 
to  do  with  difficulties  xn  reducing  the  data,  which 
are  measured  by  an  electron-beam  technique. 

3.  Good  agreement  with  measured  values  of  nN2  is  obtained 
at  downstream  stations.  However,  at  the  station 
closest  to  the  nozzle,  significant  differences  are 
observed  on  the  jet  side  of  the  mixing  layer,  where 
the  measured  N2  density  is  significantly  higher  than 
the  calculated  value.  The  reason  for  the  discrepancy 
is  not  clear  at  present.  If  this  increased  mixing 
were  accompanied  by  increased  temperature  spread,  it 
could  have  a  significant  effect  on  plume  emission  in 
this  region. 

> 

The  comparison  does  not  provide  a  conclusive  validation  of  the 
calculation,  because  of  questions  regarding  data  reduction  and 
influence  of  the  forebody  on  the  plume  nose  region.  However, 
the  semi -quantitative  agreement  obtained  indicates  that  the 
calculation  is  at  least  reasonably  correct,  with  the  excep- 
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tion  of  the  shocks  themselves  and  the  inner  portion  of  the 
mixing  zone  close  to  the  nozzle.  A  comparison  of  measured 
and  predicted  temperatures  would  be  most  desireable,  but  no 
comparison  can  be  made  at  this  time. 

2.2  Further  flow  field  calculations ; 

Since  our  last  report,  a  number  of  additional  calcula¬ 
tions  of  Atlas  and  Apollo  exhaust  plume  flow  field  calcula¬ 
tions  have  been  performed.  These  calculations  (and  those 
reported  in  the  proceeding  section)  were  made  at  the  AEC 
Lawrence  Laboratory  in  Berkeley,  California,  using  the 
Laboratory's  CDC.  7600  computer.  In  addition  to  being  of 
interest  in  their  own  right,  these  calculations  allowed  us 
to  gain  some  experience  with  the  Lawrence  Laboratory  sys¬ 
tem,  which  we  expect  to  be  using  in  most  future  work  on  this 
contract . 

Computations  performed  included  a  laminar  Atlas  sus- 
tainer  plume  at  180  km,  a  "turbulent"  oustainer  plume  at 
120  km  (using  the  same  eddy  viscosity  model  used  in  our  pre¬ 
vious  Apollo  calculation) ,  and  two  laminar  Apollo  S-II  plumes 
at  180  km,  one  at  the  minimum  velocity  at  this  altitude 
(essentially  the  orbital  injection  altitude)  of  4.2  km/sec 
and  one  at  che  burnout  velocity  of  5.3  km/sec.  Because  of 
an  input  data  error  in  the  continuation  of  the  "slow"  Apollo 
calculation,  only  the  first  3  km  of  this  plume 
calculation  is  reliable.  It  was  our  intention  to  perform  a 


number  of  calculations  of  chemistry  and  excitation  processes 
in  these  flow  fields  at  our  Detroit  facility.  However,  a 
number  of  difficulties  arose  in  transferring  the  taped  output 
from  these  calculations  to  regular  CDC  tapes,  and  although 
these  have  now  been  solved  we  have  as  yet  not  undertaken  any 
calculations  of  rate  processes  in  these  flow  fields.  We  plan 
to  do  so  in  the  next  reporting  period. 


Multistep  V-V  and  V-T  processes: 


In  the  previous  report,  we  began  calculations  of  the 
effects  of  multiple-step  excitation  processes  upon  plume 
emission.  These  processes  involve  the  production  of  poly¬ 
atomic  molecules  in  states  with  large  vibrational  quanta, 
such  as  the  v3  mode  of  HjO,  through  lower-energy  intermediate 
states,  'such  as  the  HjO  \>2  mode.  The  reason  that  such  pro¬ 
cesses  are  of  interest  is  that  computations  of  plume  flow 
fields  indicate  that  moat  of  the  HjO  in  the  exhaust  plume  is 
cold.  Available  information  upon  the  rates  of  direct  excita¬ 
tion  processes  for  the  HjO  and  v3  modes  (in  the  form  of  rate 

constants  derived  from  thermally  averaging  Marriott's  quantum- 

6 

mechanical  cross-sections)  indicates  that  these  are  slow  at 
low  temperatures.  This  behavior  is  associated  with  the  size 
of  tho  vibrational  quantum  which  must  be  produced  in  the  colli¬ 
sion.  As  a  result  of  the  plume’s  being  cool  and  the  direct  ex¬ 
citation  rate's  being  slow,  it  is  possible  that  an  indirect  ex- 
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citation  mechanism  involving  two  smaller  quanta  exchanged  in 
two  separate  collisions  may  be  faster  than  the  single  quantum, 
one-collision  process. 

The  scaling  arguments  introduced  in  our  previous  report 
indicate  that  these  multiple-collision  processes  should  be 
less  important  at  high  altitude,  since  (at  constant  tempera¬ 
ture)  the  volume -integrated  rate  should  vary  as  p1>/4  to  p3y^. 
depending  upon  the  rate  of  the  second  process.  In  addition, 
some  comparative  calculations  of  the  rates  of  a  particular 
two-step  vibrational  excitation  process  and  the  chemi-excita- 
tion  process  considered  earlier  indicated  that,  in  the  S-II 
plume  at  120  km,  there  should  not  be  substantial  differences 
in  total  rates  of  production.  The  process  considered  was  the 
V-T,  V-V  process 

1.  H20  +  M  3=£H20(V2)  +  M 

2.  H20(v2)  _ ►  H20  +  hv2  (6.3p) 

3.  H20  (v2)  +  M  HjO  (v^  or  v3) 

4.  N2  +  M  s=±  N2+  +  M 

5.  N2+  +  H20  (v2)  jp£  H20  (vx  or  v3)  +  N2 

6.  H20  (v3  or  v3)  — ►  H20  +  hv3  Qr  3  (2.7 m) 

in  which  the  V-V  transfer  in  the  4^  step  is  close  to  resonance 
and  snould  occur  frequently.  In  the  current  reporting  period, 
we  have  examined  this  process  in  greater  detail;  however,  using 
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the  rate  data  available  to  us,  we  find  no  great  contribution 
to  the  radiant  intensity  from  this  mechanism. 

Radiative  losses  from  the  H20  v2  mode  (step  2)  are  sig~ 
nificant  enough  that  this  step  introduces  some  stiffness  into 
the  system  of  differential  equations  giving  the  rates  of 
these  processes.  In  previous  applications  all  rates  were 
slow  enough  that  no  problems  of  this  kind  were  encountered, 
and  it  became  necessary  to  modify  slightly  the  numerical  in¬ 
tegration  procedure  in  the  NONEQ  code  to  treat  this  problem 

satisfactorily.  In  this  exploratory  work,  we  have  adopted 

7 

a  procedure  similar  to  that  of  Moretti,  though  much  simpler 
inasmuch  as  only  one  stiff  step  occurs. 

We  assume  that  the  radiative  depletion  of  the  and 
modes  is  fast  enough  that  we  may  neglect  the  reverse  processes 
in  steps  3  and  5.  The  rate  of  photon  production  (step  6)  is 
then  essentially  the  rate  of  production  of  H20  molecules  in 
the  or  modes,  as  all  these  molecules  radiate  immediately. 
This  assumption  causes  us  to  overestimate  the  radiant  intensity 
due  to  this  process.  Since  the  calculated  intensity  is  low, 
this  overestimate  causes  no  problem.  With  these  assumptions, 
the  rates  of  production  of  H20  (v2),  H20  (v^  or  v^) ,  and  N2+ 
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Ht  XH20(v2)  kl,f  XH20  XM  ”  (kl,r  XM  +  k2  +  k!5,f  XN2t)XH20(v2) 

(2.1) 

3t  XN2t  =  k4,f  XN2  XM  "  (k4,r  XM  *  k5,f  XH20(v2))xN2f  (2.2) 

di 

dt  XH20(v1  or  v3)  =  k3,f  XM  xH20(v2)  *  k5,f  xH20(v2)  xN2t(2*3) 

Here  terms  involving  M  are  to  be  interpreted  as  sums  over  all 
collision  partners. 

Let  xh20(v2)=  x-  xN2t“  y*  »nd  xh20(v1  or  v3)  *  z  “d 

assume  that  the  mole  fractions  of  H20,  N2 ,  and  M  are  essen¬ 
tially  unchanged  by  this  set  of  processes.  Then  the  above 


equations  become 

if  -  *-i 

B  +  Cyjx 

(2.4) 

&  -  »-| 

E  +  Cxjy 

(2.5) 

it  -  [** 

CyJ  x 

(2.6) 

where  the  coefficients  A  to  F  are  slowly  varying  function*,  of  time 
along  a  streamline.  Since  y  also  varies  slowly  (the  N2  vibra¬ 
tional  relaxation  times  are  long) ,  over  a  small  finite  time  in¬ 
terval  we  can  replace  y  by  yQ,  its  value  at  the  beginning  of 
the  interval.  Then 


-  *  -  [®  +  =y0)x 


and  the  solution  for  x  at  the  end  of  an  interval  At  is 


~  (B  +  Cy  )  At  v 

A{X  '  e  }  +x  e-(B+Cy0)4t 


X  - 


(2.7) 


B+Cy, 


When  B  is  large,  corresponding  to  significant  radiative  losses 
from  the  x>2  m°de,  the  second  term  dominates.  Now  if  y  changes 
by  only  a  small  amount  in  time  At,  substituting  the  second 
term  from  (2.7)  in  (2.5),  letting  y  85  yQ  on  the  right,  and 
integrating  between  t  *  0  and  t  *  At  gives 


Y  -  yo  =  (D  -  Eyo)At  - 


x  v  f  -  (B  +  Cy  )  At 

irzfc  [ 1  ~ e  (2*8) 


and  in  similar  fashion  we  have 


z  -  z 


F  +Cyt 
B'  +  Cv 


1  -  e 


“(B  +  Cy  )  At 


(2.9) 


When  radiative  losses  from  the  \>^  mode  are  severe,  and  produce 
considerable  stiffness  in  the  equations,  we  may  use  Equations 
(2.7),  (2.8),  and  (2.9)  to  calculate  the  changes  in  x,  y,  and 
z  over  a  single  step.  When  the  radiative  losses  are  minimal, 
there  is  no  problem  of  stiffness  and  the  straightforward  numer¬ 
ical  integration  of  Equations  (2.1),  (2.2),  and  (2.3)  can  be 
carried  out  as  before. 


We  have  applied  this  procedure  to  estimate  the  radiant 
intensity  from  the  (turbulent)  Apollo  15  plume  due  to  this 
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mechanism  at  an  altitude  of  120  km.  The  flow  field  was  com¬ 
puted  in  the  previous  reporting  period.  Hates  of  V-T  pro¬ 
cesses  involving  the  H^O  modes  were  taken  from  the  calcula¬ 
tions  of . Fisher, ^  who  has  thermally  averaged 
Marriott's  theoretical  cross-section  for  a  variety  of  pro¬ 
cesses.  Using  these  data  underestimates  the  production 
of  I^O  (v2)  in  H20  -  H20  collisions  at  the  lower  temperatures; 
however,  it  is  difficult  tc  extrapolate  the  experimental 
data,  and  the  differences  (about  a  factor  of  5)  between  the 
theoretical  and  experimental  rates  for  this  process  are  not 

of  great  significance  in  view  of  the  results  of  the  calcu- 

8 

lation.  Marriott  has  calculated  cross-sections  for  H20 
excitation  in  collisions  wi*-h  H2  and  H20  only.  These  should 
be  the  faster  processes  wherever  the  concentrations  of  H2 
and  H2O  are  high.  In  the  outer  regions  of  the  plume,  exci¬ 
tation  by  collision  with  0,  02 ,  and  N2  is  possible.  The 
rates  of  excitation  by  C>2  and  N2  should  be  slow  because  of 
the  greater  mass  of  these  molecules,  but  the  rate  of  excita¬ 
tion  by  O  -  H2O  collisions  may  be  important.  In  the  absence 
of  any  theoretical  or  experimental  information,  we  have 
assumed  (as  in  our  previous  report)  that  these  rates  are  iden¬ 
tical  to  the  H20  -  HjO  rates,  for  which  the  reduced  mass  of 
the  collision  pair  is  similar. 

The  N2  v-t  rates  are  derived  from  the  correlation  of 
9 

Millikan  and  White.  These  should  be  satisfactory  for  N2  and 
02  collisions,  but  their  use  for  H2,  H^O,  and  0  is  question- 


1 


I 


able.  No  information  is  available  on  relaxation  of  N2  by  H2 
or  H20;  however,  the  relaxation  of  the  electronically  similar 
CO  molecule  by  H2  is  part  of  the  data  on  which  the  correlation 
is  based,  and  extension  to  N2  appears  reasonable.  The  relaxa¬ 
tion  of  N2  by  0  atoms  is  currently  under  study  at  Stanford 
Research  Institute,  but  no  results  were  available  to  us  at 
the  time  these  calculations  were  being  done.  We  assume  that 
the  probability  of  transferring  a  vibrational  quantum  between 
N2  and  H20  is  0.01  at  all  temperatures;  no  data  are  available, 
but  on  the  basis  of  N2  -C02  data  this  value  would  appear  to 
be  a  somewhat  higher  estimate. 


It  is  generally  agreed  that  the  vibrational  temperature 
of  N2  is  not  in  equilibrium  with  with  the  translational  tem¬ 


perature  at  altitudes  above  100  km.  A  recent  calculation 
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by  Jamshidi  and  Kummler  (Figure  1)  gives  a  value  of  approxi¬ 


mately  7  50 °K  for  Tv  M  at  120  km.  The  translational  temper- 

V  ,  N2 


ature,  according  to  the  CIRA  tables,  is  355 °K.  In  our  cal¬ 
culations  for  the  radiant  intensity  of  the  S-II  stage,  we 


have  used  T„  »  750°K  and  „  *  2000°K.  The  higher  tern- 

V|N^  V  #N  2 

1  1  ii 

perature  corresponds  roughly  to  a  calculation  by  Walker  assum¬ 


ing  that  all  energy  deposited  in  the  1D  state  of  atomic  oxygen 
by  solar  photodissociation  of  02  (Schumann-Runge  continuum) 


is  channeled  into  Njt. 


Results  of  these  calculations  are  shown  in  Figure  2* 
Here  we  show  cumulative  radiant  intensity  as  a  function  of 
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!  1 


■  -jfia 


axial  distance  behind  the  vehicle  for  the  two  values  of 

T..  M  .  As  expected,  the  calculated  intensity  varies  strongly 

v  ,«2 

with  ambient  vibrational  temperature.  That  it  does  not  vary 
even  more  strongly  is  a  consequence  of  excitation  and  de-ex¬ 
citation  of  N~  +  in  the  shock  layer.  At  a  value  of  T,.  ..  = 

2  *  V,N_ 

** 

750 °K,  additional  N2  +  is  produced  in  the  forward  3  km  of  the 

plume  shock  layer;  however,  at  Ty  N  =  2000°K,  relaxation 

2 

processes  in  the  shock  layer  actually  decrease  the  amount  of 
N^t,  since  most  of  the  shock  layer  is  cooler  than  the  ambient 
vibrational  temperature.  The  rates  of  production  are  shown 
(again  cumulatively  as  a  function  of  downstream  distance)  in 
Figure  3  .  The  slope  of  these  curves  gives  the  local  produc¬ 
tion  rate  per  unit  length  of  the  plume  shock  layer. 

The  values  of  total  radiant  intensity  from  the  h20 

and  v3  bands  are  1.5  x 10^  watts/steradian  at  Ty  N  =  750°K 

4  ^ 

and  3.5  x  io  watts/steradian  at  Tv  M  =  2000 °K.  These  values 

v  •  n2 

would  be  up  to  a  factor  of  5  higher  if  we  increased  the  rate 
constants  of  step  1  to  bring  the  H20  (v2)  production  rate  by 
collisions  with  H20  into  agreement  with  experimental  values. 
Recall,  however,  that  most  of  the  other  assumptions  that  we 
made  tended  to  favor  the  production  of  2.7m  radiation  from 
this  system.  It  is  therefore  difficult  to  ascribe  the  S-II 
emission  to  this  mechanism  unless  one  or  more  of  the  important 
rate  constants  were  in  error  by  a  factor  of  100 . 


2 . 4  IR  chemiluminescence  in  hydrocarbon/oxygen  systems; 


In  our  previous  report,  we  presented  results  for  cal¬ 
culations  of  possible  IR  chemiluminescence  from  the  Atlas 
sustainer  and  Apollo  S-II  plumes.  These  calculations  were 
based  upon  a  reaction  sequence 

H2  +  0  ^  OH  +  H 

OH  +  R-H  — ►  R*  +  H20i 

where  R  is  anything  from  an  H  atom  to  a  hydrocarbon  chain 
and  unit  quantum  yield  is  assumed  in  the  second  step.  The 
results  of  the  calculations  were  generally  negative;  on  the 
basis  of  available  rate  data,  the  production  rate  of  excited 
H20  by  this  mechanism  was  quite  low. 

This  conclusion  was  surprising  in  view  of  earlier  ap- 
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proximate  calculations  by  Fisher,  who  had  indicated  that,  for 
quantum  yields  near  unity,  the  production  rate  should  be 
high.  As  a  consequence,  part  of  our  effort  during  this  report¬ 
ing  period  has  been  devoted  to  checking  the  numerical  integra¬ 
tion  procedure  in  various  ways  and  to  determining  why  our  n 
merical  calculations  differ  from  Fisher's  approximate  calcula¬ 
tion.  We  have  as  yet  found  no  indication  that  the  numerical 
calculations  were  incorrect.  In  general,  any  approximations 
made  tend  to  favor  the  production  of  excited  H20  rather  than 
to  limit  its  production.  Fisher's  calculation  differs  from 
the  current  calculation  in  three  ways;  1)  he  assumed  a  higher 


value  for  the  rate  constant  of  the  0  +  H0  reaction  than  that 
currently  favored  (Leeds  compilation)  ;  2)  the  extrapolation  | 

J 

which  he  introduced  to  extend  the  data  from  the  two  streamlines 
available  to  him  at  that  time  overestimates  the  total  reaction; 

I 

and  3)  he  assumed  the  concentration  of  CH,0,  which  turns  out 

z  i' 

! 

to  be  the  most  important  R-H  species  in  hydrocarbon/oxygen 

i 

plumes  because  of  the  high  value  for  the  OH  +  CH00  rate  con-  I 

z  j 

stant,  to  be  about  2%  by  moles  everywhere  in  the  combustion  j 

gas,  whereas  in  the  current  calculation  it  is  assumed  to  be  I 

J 

this  high  only  in  the  turbine  exhaust  gases. 


It  is  difficult  to  assign  mechanistic  implications  to  j 

i 

these  calculations  for  hydrocarbon/oxygen  systems  because  it  •  j 

has  been  necessary  to  guess  the  amounts  of  R-H  species  (other  j 

! 

than  H2)  which  are  present  in  the  plume.  Ther  are  no  avail- 

( 

able  measurements  of  these  concentrations.  All  available  cal- 

i 

culations  are  based  upon  assumptions  of  local  microscopic  1 

homogeneity  and  thermochemical  equilibrium  in  the  combustion 
chamber,  and  while  these  assumptions  lead  to  reasonable  values  j 

of  exit  conditions  and  major  species  concentrations,  they  are  j 

S 

) 

wholly  inadequate  for  predicting  minor  species  levels.  Another  > 

j 

serious  defect  is  that  the  reaction  system  is  incomplete.  A  i 

j 

potentially  large  number  of  important  reactions  was  purposely 

i 

excluded  from  the  calculation.  Hydrocarbon  combustion  is  a 
very  complex  process,  and  to  include  all  possible  reactions  in 
a  calculation  of  radiant  intensity  would  prove  extremely 
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lengthy  and  expensive.  (On  the  order  of  1000  different  reac- 
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tions  occur  in  the  oxidation  of  a  gasoline-fraction  alkane.  ) 
For  most  of  these  reactions,  rate  data  are  unavailable  and 
must  be  estimated;  thus  increasing  the  complexity  of  the  cal¬ 
culation  does  not  necessarily  increase  its  accuracy.  Our 
decision  to  limit  the  calculation  to  essentially  the  above 
sequences  was  based  upon  the  outcome  of  the  approximate  cal¬ 
culations,  which  suggested  that  this  sequence  was  sufficient 
to  explain  observations.  It  now  appears  that  this  assumption 
is  incorrect,  and  that  the  reaction  scheme  should  be  expanded. 

14 

Fisher  and  Kummler  have  reviewed  hydrocarbon  ^oxidation 
processes  as  part  of  studies  directed  at  the  interpretation 
of  flow-tube  experiments  in  IR  chemiluminescence.  The  major 
features  of  the  combustion  process  ir.  the  presence  of  sig¬ 
nificant  amounts  of  0  atoms  are  the  following:  A  significant 
contribution  to  the  total  OH  production  rate  is  the  attack  of 
0  atoms  upon  secondary  hydrogen  atoms  in  aliphatic  hydrocar¬ 
bons  : 

o  +  r1-ch2-r2  -  OH  +  R1“CH-R2 

This  reaction  is  generally  nearly  thermoneutral,  so  that  no 
excited  products  are  expected.  The  radical  can  react  with 
either  0  or  OH  to  form  an  olefin;  if  R2~  is  R3-CH2~  ,  then 

R1-Oi-CH2-  R3  +  0  — ►  R1-CH»CH-R3  +  OH 

r1-ch-ch2-r3  +  oh  — ►r1-ch-ch-r3  +  h2o 
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Enough  energy  is  liberated  in  these  reactions  to  excite  the 
products.  The  olefin  in  turn  can  react  with  0  to  form  an 
aliphatic  radical  and  the  formyl  radical: 

R1-CH=CH-R3  +  0  — ►  R1-CH-R3  +  CHO 

Then  the  formyl  radical  can  undergo  reactions  with  0  to 
give  CO  and  OH,  with  OH  to  give  0  and  formaldehyde,  or  with 
H20  to  give  OH  and  formaldehyde.  A  series  of  reactions  in¬ 
volving  0,  OH,  H2,  and  CO  completes  the  system. 

We  are  currently  coding  a  lumped -parameter  reaction 
scheme  for  this  sequence  of  processes.  We  expect  to  report 
on  the  results  of  these  calculations  in  our  next  semi-annual 
report.  While  the  scheme  was  being  formulated,  we  performed 
an  approximate  calculation  in  which  we  merely  added  to  the 
reaction  scheme  used  in  our  previous  Atlas  sustainer  calcula¬ 
tions  the  reaction 


n  -  C5H12  +  0 


C5H11  +  0H 


assuming  that  all  the  hydrocarbons  in  the  sustainer  turbine 
exhaust  flow  were  normal  pentane .  The  rate  constant  measured 
by  Herron  and  Huie15  is  8  x  1013  e“2320/T  CIn3ymole  aecona. 

We  attempted  to  account  for  the  fact  that  there  is  a  chain 
of  subsequent  processes  by  ignoring  the  depletion  of  the  hydro¬ 
carbon  species  in  this  reaction.  Results  of  this  calculation 
for  slightly  different  initial  formaldehyde  concentrations  in 
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the  turbine  gas  are'  shown  in  Figures  4  and  5,  together  with 
results  from  a  similar  calculation  without  this  reaction. 
Including  OH  production  by  this  process  greatly  increases  the 
potential  level  of  IR  chemiluminescence  in  the  sustainer 
plume,  and  makes  this  mechanism  again  appear  a  likely  con¬ 
tributor  to  the  overall  radiant  intensity  in  the  OH  and  H20 
vibration-rotation  bands.  More  detailed  calculations,  in¬ 
cluding  variation  of  the  assumptions  about  the  concentrations 
of  unburned  fuel  and  combustion  intermediates,  are  necessary 
before  we  can  draw  any  more  firm  conclusions. 

This  mechanism  will  not  contribute  to  emission  from 
hydrogen-oxygen  plumes  such  as  those  of  the  Apollo  S-II  or 
Centaur  vehicles.  It  may  contribute  to  emission  in  amine- 
fueled  vehicles,  although  it  seems  unlikely  that  many  long- 
chain  hydrocarbons  are  produced  in  the  combustion  of  methyl¬ 
ated  hydrazines.  There  may  be  analogous  reactions  in  sub 
stituted  hydrazine  systems;  we  are  currently  conducting  a 
literature  survey  to  discover  whether  this  is  the  case. 
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3.0  Simulation  Studies 


During  this  reporting  period  we  initiated  studies  of 
the  simulation  of  high-altitude  exhaust  plume  emission  in 
ground  facilities  and  in  high-altitude  atmospheric  releases. 
These  studies  are  aimed  at  simulating  the  plume -atmosphere 
mixing  region  rather  than  the  plume  core,  which  can  be  re¬ 
produced  (except  for  very  large  engines)  in  ground  facilities. 
There  seem  to  be  two  reasons  for  wanting  to  perform  a  simu¬ 
lation  experiment: 

1)  The  number  of  processes  which  could  contribute  to 
the  plume  emission  in  the  near  IR  is  large,  and 
the  processes  responsible  have  as  yet  not  been 
clearly  identified.  A  simulation  experiment  of¬ 
fers  (at  least  in  concept)  the  chance  to  repeat 
the  experimental  conditions,  to  ask  more  detailed 
questions  about  the  state  of  the  plume  as  a  func¬ 
tion  of  position,  and,  by  changing  the  experimen¬ 
tal  conditions,  to  alter  the  rates  of  competing 
processes  to  determine  which  are  most  significant:. 

In  particular,  it  is  desireable  to  find  out  whether 
any  important  effects  have  been  left  out  of  the 
calculation . 

2)  The  models  developed  for  predicting  plume  emis¬ 
sion  are  complicated,  and  verifying  the  computa¬ 
tional  results  is  difficult.  A  simulatior  ex¬ 
periment  offers  the  chance  of  exercising  the 
models  against  a  representation  of  reality  in  which 
more  detailed  comparisons  are  possible  than  are 
available  from  most  field  data. 
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Both  ground  and  atmospheric  simulation  experiments  have  some 
advantages  and  some  disadvantages  in  terms  of  achieving  these 
objectives. 

The  advantages  of  a  ground  simulation  experiment  are 
its  controllability,  reproducibility,  capacity  to  accept 
changes  in  experimental  conditions ,  and  the  fact  that  a  wide 
variety  of  diagnostic  measurements  (local  and  remote)  can  be 

i 

performed.  One  can  simulate  both  proper  ratios  of  binary  reac- 
tion  rates  and  plume  fluid  dynamics  at  the  same  time  if  the  product 
of  engine  thrust  and  ambient  pressure  is  held  constant  and  all  tem¬ 
peratures,  velocities,  and  pertinent  compositions  are  reproduced. 

However,  one  cannot  simulate  effects  depending  upon  the  ratio 
of  the  radiative  lifetime  to  a  characteristic  collision  time 
(or  flow  time)  in  the  infrared  if  one  scales  the  thrust  by 
more  than  about  one  order  of  magnitude.  The  simulation  can 
break  down  in  either  of  two  ways: 

1)  The  increased  pressure  can  result  in  the  radiation's 
being  forced  toward  the  equilibrium  limit,  in  which 
the  local  emission  is  independent  of  the  processes 
producing  the  excited  states  and  the  experiment 
tells  one  nothing  about  these  processes. 

2)  The  decreased  length  means  that,  even  if  radiative 
decay  is  the  predominant  de-excitation  mechanism, 
excited  molecules  can  be  convected  into  downstream 
parts  of  the  flow  field  before  they  are  de-excited. 

Since  the  radiation  is  no  longer  characteristic  of 
local  conditions,  its  interpretation  is  difficult. 


In  either  case,  the  experiment  cannot  be  scaled  back  to  the 
real  plume  conditions  without  considerable  effort.  Other 
difficulties  arise  when  one  attempts  to  include  full  chemical 
simulation  of  the  atmosphere  (degree  of  oxygen  dissociation, 
N2  vibrational  temperature)  and  the  nozzle  exit  gases  (com¬ 
bustion  intermediates) . 

In  the  ultraviolet  and  visible  the  difficulties  with 
the  radiative  lifetime  do  not  exist,  since  the  lifetimes  of 
all  the  upper  states  of  interest  are  so  short  that  tr/tc  <<  1 
in  both  the  real  and  the  simulated  flows.  The  major  problems 
are  associated  with  reproducing  the  nozzle  exit  chemistry. 

In  the  current  reporting  period  our  activity  has  been 
mostly  concerned  with  defining  the  possibilities  of  partial 
simulation  of  exhaust  plumes  by  means  of  atmospheric  releases 
of  material.  The  release  could  take  several  forms,  but  the 
end  result  is  essentially  a  cloud* of  material  which  inter¬ 
acts  with  the  ambient  atmosphere  in  a  manner  which  under 
some  conditions  closely  mimics  the  interaction  of  the  plume 
with  the  atmosphere.  The  morphology  of  the  two  flows  (in 
the  continuum  regime)  is  shown  in  Figure  6  . 

We  assume  that  the  release  is  either  cylindrically  or 
spherically  symmetric,  that  its  initial  velocity  is  zero, 
and  that  the  container  is  removed  at  t»0.  The  released  gas 
expands  in  a  manner  determined  by  the  initial  distribution  of 
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density  and  the  value  of  y  behind  a  front  which,  in  vacuo, 
would  attain  a  speed  of  2aQ/(y  -  1) ,  where  aQ  is  the  initial 
speed.16  (The  density  and  pressure  at  this  front  are  zero.) 

This  front  is  the  analog  of  the  limiting  line  for  the  Prandtl- 
Meyer  expansion  of  the  plume  gas.  As  the  released  gas  ex¬ 
pands,  its  dynamic  pressure  decreases.  When  the  dynamic 
pressure  of  the  released  gas  becomes  comparable  to  the  am¬ 
bient  pressure,  the  released  gas  will  be  decelerated  by  the 
surrounding  atmosphere.  Since  most  of  the  expanding  gas 
will  become  supersonic  if  the  expansion  ratio  is  large,  the 
decelaration  will  occur  through  a  shock  wave  which  propa¬ 
gates  inward  from  the  contact  surface  into  the  expanding  gas 
cloud.  This  shock  wave  will  eventually  reach  the  center  or 
axis  of  the  cloud  and  be  reflected  (regularly  -  there  is  no 
analog  of  the  Mach  disc) .  Since  the  released  gas  will,  if 
initially  hot,  expand  at  velocities  supersonic  with  respect 
to  the  ambient  atmosphere,  the  ambient  gas  just  outside  the 
contact  surface  will  be  swept  outward  at  a  supersonic  velo¬ 
city  and  a  shock  wave  will  form  in  the  ambient  gas  (just  as 
in  a  shock  tube).  In  general,  the  contact  surface  will  over¬ 
shoot  its  equilibrium  position  and  oscillate  around  this 
position  until  the  internal  shock  reflects  several  times  off 
the  center  or  axis  and  the  contact  surface.  The  equilibrium 
radius  is  proportional  to  (E0/P„) ,  where  EQ  is  the  initial 
energy  (per  unit  length  in  the  cylindrical  case)  of  the  re¬ 
leased  material,  PM  is  ambient  pressure,  and  j  is  a  metric 
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coefficient  (1  for  cylindrical  and  2  for  spherical  geometry) . 
Viscous  transport  occurs  across  the  interface,  which  can  be 
Rayleigh-Taylor  unstable  if  the  released  gas  at  the  interface 
is  more  dense  than  the  shocked  air  (normally  the  case,  since 
the  released  gas  is  relatively  cool) . 


In  the  exhaust  plume  (cylindrical  geometry)  a  roughly 
similar  series  of  processes  takes  place.  There  is,  in  fact, 
a  direct  analogy  between  flows  in  an  exhaust  plume  with  high 
exit  and  free-stream  Mach  numbers  and  a  cylindrical  release. 
The  analogy  becomes  exact  over  the  outward -moving  portions 
of  the  flows  (forward  part  of  the  plume  and  early  part  of 
the  release)  in  the  limiting  case  that  all  Mach  numbers  are 
high,  all  flow  angles  are  small  (but  M  sin$  <<  1),  and  the 
jet  and  free-stream  limiting  velocities  are  identical.  In 
this  case  the  two  flows  may  be  connected  through  a  simple 
Galilean  transformation  x  »  uwt.  (The  aerodynamicist  calls 
this  the  "hypersonic  equivalence  principle.") 


7 


In  practical  situations,  the  conditions  for  exact  analogy 
are  not  fulfilled.  As  we  noted  parenthetically  above,  the 
plume  flow  will  form  a  Mach  disc  and  the  release  flow  cannot. 

If  the  Mach  number  at  the  nozzle  exit  is  only  moderately  large, 
the  Prandtl-Meyer  turning  angle  is  usually  large  enough  that 
the  air  shock  is  detached  at  the  plume  nose,  and  in  all  prac¬ 
tical  cases  the  small-angle  criterion  is  not  fulfilled.  In 
addition,  a  point  release  would  usually  be  spherical  (or  some 


approximation  thereto)  and  the  exact  analogy  breaks  down  in 
any  case  because  of  different  divergence  terms  in  the  equa¬ 
tions  of  motion.  Nonetheless,  there  may  remain  a  fair  degree 
of  similarity  between  the  two  flows,  and  in  particular  for 
cases  when  vehicle  velocity  and  plume  velocity  are  closely 
matched  the  development  of  a  stationary  release  as  a  func¬ 
tion  of  time  provides  a  semi-quantitative  analog  to  the  de¬ 
velopment  of  a  plume  as  a  function  of  downstream  distance. 

The  quantity  of  material  to  be  released  is  a  question 
of  considerable  importance,  since  it  sizes  the  carrier  ve¬ 
hicle  required.  If  we  take  the  approach  that  we  maintain 
constant  altitude  (to  get  identical  atmospheric  conditions) 
and  reproduce  the  plume  Knudsen  number  based  upon  its  equi¬ 
librium  radius  with  the  similar  Knudsen  number  of  a  sta¬ 
tionary  spherical  release  (to  reproduce  mixing  phenomena  or 
molecular  distribution  functions) ,  then  the  amount  of  mater¬ 
ial  can  be  determined  approximately  as  follows:  The  equi¬ 
librium  radius  of  the  plume  is 

rp  -  0.1  yk/P„  (3.1) 

and  of  the  release  is 

rR  ~  <  VP->  1/3  '  '  (3'2) 

We  assume  that  the  velocity  of  the  plume-generating  vehicle 
is  equal  to  the  limiting  velocity  of  the  plume  gases  up.  Ther. 


the  temperatures  in  the  plume  are,  for  the  most  part,  governed 

by  the  nozzle  exit  temperature  (which  provides  the  energy  which 

will  be  converted  into  kinetic  energy  of  lateral  expansion)  and 

the  air  temperature.  The  rocket  thrust  is 

P*ffl  Up  (3.3) 

and  the  initial  energy  of  the  released  material  is 

E  =  27  “  UR  (3*4) 

where  uR  is  the  limiting  steady-flow  velocity  of  the  released 
gas.  Then  if  is  the  same  in  both  cases  and  rp/rR«l,  we 
have 

0.1(m  up/PJ1/2  «  (b  uJ/2yPJ1/3 
or 

m  «  0.002y(m3  u^/P^u4) 1/2  . 

For  a  rocket  with  105  lb  thrust  at  an  altitude  of  180  km, 

o  5 

in  *  1.5  x  10  gm 

5 

up  w  3.5  X  10  cm/sec 

-3  2 

Pa  at  2  x  io  dynes/cm 

Assuming  that  the  released  gas  is  initially  at  the  nozzle 
exit  temperature,  and  taking  this  temperature  as  1/4  the  cham- 

C 

ber  temperature,  gives  uR  m  1/75  *  10  cm/sec.  Then  with 
Y  ■  1.4,  we  have 


m  «  2.8  x  10“3(3.38  xlO15  x4.3  x  1016/2  x  10  3  x  9.4  xl020)1/2 


or 

4 

m  *  2.5  x  10  gm  m  55  lb 

This  appears  to  be  within  the  lifting  capability  of  available 
sounding  rockets,  even  if  a  sizeable  extra  factor  is  included 
for  storage  and  controls. 

It  is  not  necessary  to  perform  the  release  at  zero 
velocity.  By  releasing  the  material  with  appreciable  velo¬ 
city,  one  obtains  a  higher  collision  energy  between  release 
and  atmospheric  molecules.  However,  one  loses  both  symmetry 
(since  forward  and  rear  portions  of  the  released  gas  feel  a 
different  effective  pressure)  and  the  qualitative  analogy  to 
the  development  of  a  thrusting  plume .  The  forward  part  of 
the  release  qualitatively  resembles  a  retro  plume  during  a 
portion  of  the  expansion.  Velocities  available  depend  upon 
the  choice  of  carrier  vehicle  and  range  from  essentially 
zero  (sounding  rockets)  to  about  7  km/sec  (piggyback  on 
large  missile) . 

Since  the  purpose  of  performing  the  release  is  to  learn 
about  chemistry  and  excitation  processes  in  high-altitude 
plumes,  one  hopes  that  thu  distribution  of  collision  energies 
in  a  release  can  be  tailored  to  correspond  closely  to  the 
distribution  in  a  plume.  In  this  report  we  have  begun  in- 


vestigating  the  conditions  under  which  this  correspondence 
exists.  We  intend  to  examine  both  high-altitude  plumes , 
where  the  flows  are  in  a  nearly  free-molecular  condition 
(after  the  initial  portion  of  the  expansion) ,  and  moderately 
low-altitude  plumes  (100-150  km) ,  where  the  flows  are  in  the 
continuum  regime.  In  the  low-altitude  situation,  we  will 
employ  a  new  one-dimensional  unsteady  flow  code  developed 
over  the  past  year  at  Wayne  State  University  for  the  release 
and  compare  the  results  to  MULTITUBE  calculations  of  the 
plume  flow.  In  the  high-altitude  situation,  we  adopt  a 
first  collision  approach,  and  the  remainder  of  this  section 
is  concerned  with  the  formulation  of  this  comparison  at 
high  altitudes. 

We  assume  that  the  rocket  engine  generating  the  exhaust 
plume  and  the  mass  of  released  material  are  large  enough  that 
the  limiting  density  distribution  in  each  flow  is  achieved  as 
a  result  of  essentially  inviscid  expansions.  The  ambient 
mean  free  path  X M  is  taken  to  be  very  long,  so  that  the  flow 
patterns  established  in  the  expansions  persist  to  largo  dis¬ 
tances  or  times.  At  some  distance  r*  or  time  t*  the  expanding 
plume  or  cloud  becomes  itself  so  rarefied  that  it  is  highly 
permeable  by  atmospheric  molecules.  (We  assume  that  Aw  is 
considerably  greater  than  the  dimensions  of  the  plume  or 
cloud  when  this  occurs.)  Under  these  conditions  the  initial 
collisions  between  injected  and  atmospheric  molecules  can 


be  treated  as  irare  and  independent  events,  and  the  velocity 
and  density  distributions  of  the  bulk  of  the  injected  mater¬ 
ial  are  unaffected  by  the  presence  of  ambient  molecules. 

Thus  these  distributions  are  those  appropriate  to  flow  into 
a  vacuum.  It  is  these  "first  collisions"  which  are  of 
especial  interest  in  problems  where  the  reaction  or  excita¬ 
tion  cross-section  is  a  strongly  increasing  function  of  en¬ 
ergy  . 


For  an  exhaust  plume  expanding  into  a  vacuum,  the  asymp¬ 
totic  density  distribution,  in  spherical  polar  coordinates 
(r,Q,<p)  centered  at  the  nozzle  (which  is  very  small  compared 
with  the  dimensions  r*  of  the  interaction  volume)  and  with 
the  z-axis  or  north  pole  pointing  along  the  thrust  vector 
(Figure  7) ,  is  approximately 


where  0^  is  the  limiting  angle  for  a  Prandtl-Meyer  expansion 
at  the  nozzle  lip.  The  effects  of  nozzle  wall  boundary  layer 


will  be  neglected.  Here  also  uB  is  the  limiting  velocity 

for  a  steady  flow,  uw  -  /Sh"  where  H  is  the  total  enthalpy. 

*►  A 

The  velocity  distribution  is  u(0)  =  u^rfO)  for  0  <  6a  . 


For  an  initially  uniform  gas  cloud  expanding  into  a 

vacuum,  the  asymptotic  density  distribution  can  be  taken 
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from  an  approximate  formulation  given  by  Hubbard: 


P  *  D(y  ,  j ) 


nlB(Y.j) 


(3.7) 


where  R  is  radial  distance  from  the  center  of  the  cloud,  t  is  time 
j  is  a  metric  coefficient  (equal  to  0,  1,  or  2  as  the  geometry  is 
planar,  cylindrical,  or  spherical),  m  is  the  mass  released  (per 
unit  area  or  length  in  the  case  of  a  planar  or  cylindrical  flow) , 


aQ  =  / ( y — 1 ) H  is  the  initial  sound  speed,  and 


u  =2  VH/(y-l) 
max  '  ' 


(3.8) 


is  the  speed  of  the  cloud  boundary.  The  coefficient:;  D  and  B 
are  given  by  the  mass  conservation  requirement 


T  fe*r) 


r  ^3  +  2  +  3  ~  -  j 

r (b  +  i)r(i  +  3.+ 


and  the  energy  conservation  requirement 


/  \i+2  t(b  2  + 

D  ■  (V)  [yjr^r]  - 7 — -7TT7\ 

L  J  r  (b  +  i)r/i  +  — -J 


(3.9) 


(3.10) 


The  velocity  distribution  is  u(R,t)  =  R/t  for  |r|  <  u  t  . 

max 

Equations  (3.9)  and  (3.10)  must  be  supplemented  by  an 
additional  equation  in  order  to  determine  the  values  of  all 
three  quantities  n,  B,  and  D.  In  the  planar  case,  an  exact 
solution  to  the  expansion  problem  is  available,  and  in  this 
solution  the  asymptotic  results  give 

n  =  2 


In  the  cylindrical  and  spherical  cases  no  such  solution  exists, 
and  one  must  look  to  numerical  solutions  for  assistance.  Hub¬ 
bard  has  chosen  to  fit  the  value  of  n  to  numerical  calcula¬ 


tions  for  the  cylindrical  case  and  to  determine  the  values  of 
B  and  D  by  Equations  (3.9)  and  (3.10).  The  value  of  n  chosen 


The  maximum  velocity  is  considerably  larger  than  the 
limiting  velocity  u,  in  a  steady  expansion#  by  the  ratio 


(3.11) 
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In  the  unsteady  expansion,  work  is  done  on  the  gas  at  the 

I 

!  leading  edge  of  the  expansion  by  the  gas  behind  it,  so 

i  that  a  transfer  of  energy  occurs  from  gas  toward  the  center 

; 

i  of  the  cloud  to  gas  at  its  edge.  In  the  steady  expansion, 

i 

i  there  is  no  leading  edge,  so  that  each  element  of  gas  tnot 

i 

i 

i  only  has  work  done  on  it  by  gas  upstream  from  it,  but  also 

does  work  on  downstream  gas.  The  average  velocity  of  the 
expanding  cloud  is 


uav  "  v/2E  -  v/2h7y  ■  \/l/y 


(3.12) 


since  energy  is  conserved. 

The  interaction  of  the  atmosphere  with  the  plume  gases 

in  best  described  in  a  coordinate  system  moving  with  the  ve- 

*+■ 

hide.  If  the  vehicle  velocity  V  is  large  compared  with  am¬ 
bient  thermal  velocities,  ambient  molecules  impinge  upon  the 
plume  along  essentially  parallel  trajectories  with  velocities 
-V.  The  highly  expanded  plume  gaEes  are  also  hypersonic,  so 
that  a  first  collision  takes  place  with  a  relative  velocity 

vr  *  u#r  -  V  ,  (3.13) 


where  r  is  a  unit  radial  vector  and  the  origin  of  coordinates 
is  at  the  engine . 

If  the  molecules  &xae  hard  spheres,  then  the  fraction  of 
molecules  penetrating  the  plume  to  a  depth  S  along  any  indi¬ 
vidual  trajectory  is 
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(3.14) 


For  a  thrusting  vehicle  at  zero  angle  of  attack,  we  consider 
that  the  trajectories  are  lines  parallel  to  the  axis  and  con¬ 
sider  penetration  to  a  distance  z  behind  the  axis.  The  in¬ 
tegral  in  the  exponent  then  becomes 


nov 

,  /  r 


M  I a  u  V Z  tan  8  J 

gp  oo  q 


v© 

f  d8'f(e/)  |/v2  +  U l~ 


2Vu  cos  e ' 
00 


(3.15) 


for  any  angle  8  =  tan  (rQ/z)  where  rQ  is  the  distance  from  the 
axis  of  the  trajectory.  For  a  retrofiring  vehicle,  an  iden¬ 
tical  expression  results  except  that  the  integral  is  taken  from 
zero  to  8. 


The  number  of  first  collisions  per  unit  volume  at  a  point 
in  the  plume  given  by  particular  values  of  z  and  0  is 


dN 

3V 


noo  T  Vc,p 


n  t  nov,. 
00  r 


(3.16) 


where  n^  is  the  ambient  density  and  vc  is  the  frequency  of 
collision  of  an  ambient  molecule  with  plume  molecules.  The 

total  number  of  such  collisions  between  z  and  z  +dz  is 

O  /max  0 

•j—  dz  -  2tt  /  dr  r  dz  «  27rz2/  de  tanf  37  dz 

{  ^  cos  8 

where  6max  is  the  lesser  of  8^  and  n/2. 

The  collision  energy  depends  only  on  8 |e  »  ^  Mv^(0)j  .  The 
rate  at  which  plume  molecules  between  z  and  z  +  dz  undergo  first 


the  release  connected  with  such  a  plume  would  most  reasonably 
be  performed  with  an  appreciable  velocity.  We  again  adopt  a 
a  cylindrical  coordinate  system  with  the  origin  at  the  center 
of  the  release  and  the  z-axis  pointing  in  the  direction  par¬ 
allel  to  the  flux  of  incoming  molecules  (Figure  8) .  We  con¬ 
sider  times  short  with  respect  to  Xu _ v  but  long  with  respect 

to  the  time  required  to  establish  the  asymptotic  density  dis¬ 
tribution  given  by  Equation  (3.7). 

Consider  an  elemental  volume  consisting  of  an  annulus 
of  thickness  dr  and  width  dz  at  a  distance  r  from  the  z-axis. 
Let  n(r,z,t)  be  the  local  density  of  uncollided  atmospheric 
molecules  at  time  t,  nR(r,z,t)  the  local  density  of  released 
molecules,  vm  the  velocity  of  the  released  gas  with  respect 
to  the  atmosphere,  vr(r,z,t)  the  relative  velocity  between 
released  and  atmospheric  molecules,  and  o  the  collision  cross- 
section.  The  rate  of  accumulation  of  uncollided  atmospheric 
molecules  within  this  elemental  volume  is 

2irr  dr  dz 

The  rate  flow  of  uncollided  atmospheric  molecules  into  this 
volume  is 

n(z,r,t)vw  2nr  dr 

and  the  rate  of  flow  out  of  the  volume  is 
n  (z  +  dz,r,t)  v,,,  2wr  dr  . 
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The  rate  at  which  coilisions  occur  is 


n(z,r,t)nDa  v„  2irr  dr  dz  . 
k  r 


By  setting  accumulation  equal  to  flow  in  less  flow  out  less 
collision  rate,  and  using  the  usual  Taylor  series  expansion 
for  n(z+dz),  we  arrive  at  the  following  partial  differential 


equation  for  n(z,r,t)t 

I?  +  Iz"  v»  *  vr(r,z,t)nR(r,z,t) 


(3.18) 


Here  the  relative  velocity  is 


v  = 

r  L 


v2  +  (r2  +  z2)/t2  -  2voo  )jx2  +  z2/t  cos  9 


(3.19) 


with 


jr2  +  2  2 


and  the  local  number  density  of  released-gas  molecules  iB 


„R  .  Dl  -^yl' 

M<v>3  2  v*L7/ 


(3.20) 


The  complicated  form  of  the  third  term  makes  a  general 
solution  of  Equation  (3.18)  difficult  to  obtain.  Two  special' 
case  solutions  can  readily  be  found,  first,  if  vs  is  zero, 
the  equation  reduces  to  an  ordinary  differential  equation  for 
n (t)  at  given  r  and  z,  with  solution 


n(t;r,z)  »  exp  {-/  nR°  (V^"  z2/tjdtj 
where  R  *  ^r2+  z^  .  We  can  write 


naa€(R,t)  (3.20) 
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e  =  exp 


D  m  N*v  rt 
o 


f  dt ' 1 ' 


-<2+j) 


max 


1  r . 

\n 

i 

vw 

(3.21) 


and  in  the  case  of  a  spherical  release  (j  *  2)  the  integral 
becomes 

3  r  r(|)r(B  +  l) 


±p*y 


7  /I  \ 

TIT"  "  Ml  '  B  +1) 


r(B  +  D 


(3.22) 


where  3„(a,b)  is  the  incomplete  Beta  Function  and  x  =  (r/u  „  t) 
*  max 

With  n  thus  expressed  as  a  function  of  space  and  time,  we  can 
proceed  to  evaluate  the  collision  energy  distribution  as  fol¬ 
lows: 


At  any  point,  the  rate  of  first  collision  per  unit 
volume  is 


<jv  *  e(R,t)n—  n^(R,t)o  ^ 


and  the  entire  rate  integrated  over  the  volume  of  released 
material  is 


N(t) 


4tt 


?a,°/ 

*  A 


,umaxt 


o 

n_a  m  N 


dR  R  e(R,t)nR(R,t)/t 


=  4 


=~T 


M  a 


uraaxt 


dR  Re (R 


^f1 -  (s-mI  o. 

\  max  / 
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1  2  2 

The  collision  energy  is  ^  yR  /t  .  The  rate  at  which  released 
and  atmospheric  molecules  undergo  first  collisions  at  time  t 
with  energies  between  E  and  E  +  dE  is 


dE 


where 


(§f)t  - 


}iR/t‘ 


and  for  a  spherical  release 


(*), 


n  a  m  N 


=  4n 


av 


M  «o  *4 


RJe(R,t)  D 


-(*)■ 


B 


and  thus 


(tt)t 


n»°  m  Nav  2 
4ff  2. - -y~2y  R2e(R,t)  D 


Mua"  t* 


•  (*T 


B 


(3.24) 


The  second  situation  in  which  an  easy  solution  is  found 


is  at  times  long  enough  that  nD  is  very  small.  Then  the 


third  term  in  Equation  (3.18)  vanishes  and  the  solution  for 


n  becomes  n  £  n„  everywhere.  Then  the  local  first  collision 

"*  <w  ■v® 


frequency  is  simply 


O 

dN 

3V  -  n_  n„  a  v. 


“R  ’r 


(3.25) 


with 


vr (r,z,t) 


fv2  ♦  iii, 
[  ■  t- 


■j —  “  2v»  £ 


I 


1/2 


<r2+z2>  <  umaxt 


(3.26) 


and  the  volume -in teg rated  rate  is 


N 


4* 


,w  y  w ,2 

nma  I  r  dr/  dz  nR(r,z,t) 


vr(r,z,t)  (3.27) 


o  o 

The  collision  energy  is 


E 


1  |.u2  1  ,, 
7  pvr  *  7  M 


..2  ,  r2  +  z2  z 

v-  +  — '  2v-  t 


(3.28) 
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or,  alternatively, 


“  7  y[v~  +-pr  "  2v~  F  cos  e]< 


(3.29) 


Likewise,  we  can  rewrite  N  in  terms  of  R  and  6  as 

/umaxt 

dR  R2/  de  sin  6  nR(R,t)vr  (R,6,t)  . 
o  o 

(3.30) 

We  can  now  rewrite  both  (3.29)  and  (3.30)  to  remove  the  ex¬ 
plicit  time  dependence,  using  the  substitution  n  =  R/u  t 

max 


and 


o 

N 


1  f  2  2  2  1 

—  s*  uv  +  u  n  -  2v  u  n  cos  0 1 

2  [  00  max  1  “  “max  1  J 


(3.31 


2rrD 


,tt 


f  d0  sin  0  [v2  +  u; 


,2  rj2  “  2v  u 
max  '  •  max 


-,1/2 

n  cos  0  (3 . 32) 


where  we  have  used  (3.20).  This  development  shows  that  the 
total  number  of  collisions  is  constant  in  time,  since  it  de¬ 
pends  upon  a  density  which  varies  as  t”3  and  a  volume  which 


varies  at  t  . 


In  this  case  the  collision  energy  depends  upon  both  n  and  0 
(or  r  and  z),  so  that  determining  dN/dE  becomes  more  complicated 
The  surfaces  of  constant  E  are  ellipsoids  of  revolution  about 
the  z-axis  (Figure  9) .  These  ellipsoids  either  will  or  will  not 
intersect  the  sphere  n  *  1  depending  upon  the  values  of  E, 
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v#/  and  uma3£.  The  values  of  E  of  interest  range  between 


iin[o,  \  U(v.-Umai.)2j 


1  2 
i-  u  (v  +  u  ) 

2  H  *>  max' 


In  general,  we  believe  it  likely  that  in  practical  situations 


v  >  u  .  so  that  E  is  not  zero  within  the  sphere  rj  =  1  • 
*  max 


The  number  of  collisions  with  energies  between  E  and 


E  +  dE  is  the  number  which  occurs  in  a  shell  bounded  by  the 


two  ellipsoids  at  energies  E  and  E+dE  and  the  surface  of 


the  sphere  n  =1.  At  any  point  on  this  surface,  the  shell 


thickness  is 


dE/ | VE 


and  the  number  of  collisions  associated  with  an  elemental 


area  dA  of  this  surface  is 


ds  01 


dN  dE 


(3.31) 
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uw  -m  '•  v'oHvsveft*  the  required  energy  then 


xjcom^ «  vh-.  vu*egrai  $*-  4A  over  the  surface  E  “ 
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To  evaluate  this  integral,  it  will  prove  easier  to  work 
in  cylindrical  coordinates.  We  introduce  the  variables 
r  =  r/u  t  and  z  =  z/u  t  so  that 


max 


E 


max 

J  m|v£  +  ui„(r2+  *2)  “  2vmu - z 


max  ~ 


*  max  - 


(3.35) 


and 


0 

dN 

<3v 


n^a  D  --—-3-  |l  -  (r2+z2)n/2J 


B 


v2,  u2  (r2  +  z2)  -  2v  u  z 
®  max':  i  '  ®  max  ~ 


1/2 


(3.36) 


The  integral  over  the  surface  &  »  constant  can  be  reduced  to 
an  integration  in  the  plane  z  ■  0  by  projecting  the  surface  on 
the  plane.  Thus  if  dA  is  the  element  of  area  on  the  surface 
E  =  constant  and  dA'  is  the  element  of  area  in  the  plane  z-  0, 
these  elements  are  related  as 

dA  ■  (0§L + ($L + T  “*  (3-37' 

2  2  2 

where  x  +  y  *  r  and  dA'  «  dx  dy  *  2irr  dr.  Since 


IRi-  f 


!l. 


:;V?5 


where  the  minus  sign  is  chosen  in  the  explicit  expression 

2  2 

for  z  because  of  the  requirement  that  z  <  1  -  r  and  the 

«■»  *W  *w 

assumption  that  vB  >  umax  .  These  expressions,  when  inserted 

into  Equation  (3.40),  give  an  expression  which  depends  only 

upon  E  and  the  other  parameters  of  the  calculation.  Since 

if  vTO  >  umax  the  tangent  planes  to  the  ellipsoids  of  constant 

E  which  lie  parallel  to  the  z-axis  touch  the  ellipsoids  at 

points  outside  the  sphere  n  =1/  the  value  of  r  is  deter- 

max 

mined  by  the  intersection  of  the  surface  of  constant  E  with 
the  sphere  n  *  1: 


Again  the  minus  sign  is  chosen  to  keep  z  and  r  less  than  1. 

For  the  general  case,  it  appears  necessary  to  solve  Equa¬ 
tion  (3.18)  by  numerical  methods,  substitute  the  resulting 
values  of  n<B(r,z,t)  into  Equation  (3.36),  and  continue  the 
subsequent  integrations  (with  na  under  the  integral  sign)  in 
order  to  arrive  at  a  collision  energy  distribution.  We  have 
set  this  problem  aside  at  the  moment,  and  are  currently  check¬ 
ing  out  a  small  computer  code  which  calculates  the  collision 
energy  distributions  in  various  situations.  We  also  intend 
to  perform  a  few  spherically  symmetric  vacuum  expansion  calcu- 


lations  using  the  CELLS  code  (a  one-d.imensional ,  time-dependent 
Lagrangian  flow  code  developed  at  Wayne  State  University)  in 
order  to  examine  Hubbard's  choice  of  n  under  these  conditions. 

These  calculations  will  give  only  a  first  approximation 
to  the  collision  energy  distribution  at  high  altitude,  since 
secondary  collisions  and  the  effects  of  ambient  temperature 
have  been  omitted.  Since  including  these  effects  requires 
a  solution  to  the  Boltzmann  equation,  by  Monte  Carlo  or  other 
procedures,  we  have  no  firm  plans  to  include  them.  (The  1-D 
Monte  Carlo  code  developed  at  Wayne  State  is  suitable  for 
the  release  geometry.)  Since  currently  available  informa¬ 
tion  on  the  cross-sections  of  the  most -import ant  V-T  excita¬ 
tion  processes  indicates  that  they  have  a  strong  dependence 
on  collision  energy,  the  first  collision  treatment  appears 
to  be  a  reasonable  first  approximation . 

In  the  lower-altitude  case,  where  continuum  equations 
of  motion  are  appropriate,  we  may  compare  particle  histories 
as  calculated  by  the  CELLS  and  MULTITUBE  codes.  These  com¬ 
parisons  will  be  forthcoming  in  the  near  future. 
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Figure  3,  Cumulative  Production  of  Vibrationally- 
Excited  No  in  the  "Turbulent"  Apollo  S-II 
Plume  at  120  km  Upstream  of  a  Given  Axial 
Distance 


Figure  4. 


x ,  km. 

Cumulative  Radiant  Intensity  of  the  Laminar 
Atlas  Sustainer  Plume  at  120  km  in  the  m 
vi  or  vj  Bands  Upstream  of  a  Given  Axial 
Distance  for  the  OH  Chemiluminescence 
Mechanism.  For  curve  I,  OH  is  formed  only 
BY  G  t  ^n^2n*»*2  At*so  contributes,  Turbine 
exhaust  composition  by  moles  is  262  Ho, 

11X  CO,  n  C02,  1.3X  H2CO,  2.6X  CHii,  26X 
C2H6,  19, 5Z  CnH2n+2 ■ 


CUMULATIVE  RADIANT  INTENSITY,  kw/ster 


Figure  7 ♦  Coordinate  System  for  Describing  Rocket 
Exhaust  Plume  Flow  and  Interactions  with 
the  Free  Stream 


Figure  8.  Coordinate  System  for  Describing  Spherical 
Releases  and  their  Interactions  with  Ambient 
Gas 
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Figure  9. 
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Surfaces  of  Constant  First  Collision  Energy 
between  Atmospheric  Molecules  and  Released 
Gases  when  the  Initial  Velocity  of  the 
Released  Material  is  Large 
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